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Evidence for the validity of the Fluctuation Theorem (FT) in an atmospheric Global Circulation 
Model is found. The model is hydrostatic with variable numbers of vertical levels and different 
horizontal resolutions. For finite time intervals the largest local Lyapunov exponent (LLLE) is 
found to be negative consistent with predictions of the FT. The effect is present for resolutions up 
to wave numbers £ = 42 (~ 250km) and 10 levels. 



PACS numbers: 92.60.Bh, 92.70.Np, 05.70.Ln 
I. INTRODUCTION 



The fluctuation theorem (FT) Q, [|| describes the 
probability of negative entropy production in nonlinear 
nonequilibrium systems. This deviation from the sec- 
ond law is found on finite time scales for small (or meso- 
scopic) systems and vanishes in the thermodynamic limit. 
Gallavotti and Cohen provided a proof of the FT for time- 
reversible Anosov dynamical systems. Dewar derived the 
FT based on a maximum entropy production principle 3[ . 
For reviews on different formulations of the FT see [1, 5| . 
The validity of the FT has been proven in systems with 
low numbers of degrees of freedom, see for example @- 
In the atmosphere violation of the second law can be as- 
sociated with an increase of predictability 0, Q . The FT 
has previously been observed to hold in experiments and 
numerical simulations of turbulent flows j^, [l(| 

The modelling of continuous systems is restricted to a 
finite number of degrees of freedom which has to be mod- 
erate if long time spans have to be covered. Therefore, 
we can expect that deviations from the second law as 
described by the FT may occur also in numerical simu- 
lations. The transient increase of predictability (denoted 
as return of skill) has been studied in low dimensional 
systems (TTI | . 

In this paper we consider simulations with an atmo- 
spheric Global Circulation Model (GCM). This type of 
models are the hydrostatic dynamical cores of complex 
weather and climate models. The model is driven by a 
temperature relaxation and is subject to linear surface 
friction and hyper-diffusion. The entropy production is 
assessed by the largest local Lyapunov exponent (LLLE) 
in close nonlinear model trajectories 12]. By changing 
horizontal and vertical resolutions the number of the de- 
grees of freedom is modified. In Sections [XT] and IIIII the 
model and the analysis are described. Section ITVl includes 
the results for the statistics of the LLLE and the validity 
of the FT, and in Section [V] we summarize the results 
and discuss the conclusions. 



II. GLOBAL CIRCULATION MODEL 

In this paper we use the model PUMA (Portable Uni- 
versity Model of the Atmosphere [l3|, (HI), a hydrostatic 
global atmospheric model based on the multi-layer prim- 
itive equations expanded in spherical harmonics. Model 
variables are vorticity, horizontal divergence, tempera- 
ture and logarithmic surface pressure. 

The horizontal resolution is varied between total spher- 
ical wave number £ = 15 and 42 (7.5° and 2.8°). The 
vertical model resolution is between L — 5 and 20 verti- 
cal levels with equidistant a = p/p s coordinate where p 
is pressure (p s at the surface). For the integration of the 
nonlinear terms all fields are transformed to grid point 
space on a Gaussian latitude/longitude grid. The wave 
numbers are restricted according to the triangular trun- 
cation which is denoted as 'TT with the total wave num- 
ber £. 

To maintain a stationary state the model is driven to- 
wards a constant temperature profile with an equator-to- 
pole gradient. Dissipation is given by Rayleigh friction in 
the boundary layer and hyperdiffusion (cx £ & ) to account 
for subscale processes and numerical stability. Orography 
is not specified. The model is integrated by a leap-frog 
method with a time step of 30min (15min for £ — 42 
resolution) . 

Due to the neglect of subscale parameterizations the 
error growth in a dynamical core like PUMA corresponds 
to the long time growth regime in complex weather and 
climate models which incorporate a rapid error growth 
caused by small and fast processes like convection [15j . 



III. ANALYSIS 

As a measure for predictability in the model the lin- 
ear stability is determined by the simulation of two close 
nonlinear trajectories. One of the simulations is denoted 
as reference trajectory x(t) while the second is consid- 
ered as a perturbed trajectory x(t) initialized with weak 
random deviations in the surface pressure. The distance 
between the trajectories ||x(t) — x(t)\\ is determined using 
an Euclidean metric for all model variables in grid point 
space. Linear instability of the reference trajectory is as- 
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sessed by maintaining small distances through a regular 
rescaling of the perturbation after a time span t = 10 
days to the initial distance do (these intervals are denoted 
as rescaling intervals in the following). To collapse the 
perturbation onto the most unstable direction the anal- 
ysis neglects an initial spin up of 10 rescaling intervals in 
the 200 years trajectories. 




LLLE [1/day] 

FIG. 1. Relative frequencies of the LLLE for r = 100 h with 
a Gaussian fit (black) for the resolution T21L5. 

Within the rescaling intervals of duration tq the growth 
of the distance d is measured every 5 hours to estimate 
the local largest Lyapunov exponent 

LLLE(t, r) = - log (1) 
t d 

Here d(t, r) = ||x(f + r) — x(t + r)||, where t = nro mea- 
sures the absolute time and r < to is the time elapsed 
after the last rescaling to the initial distance do. The 
unit for the LLLE is 1/day. In the following the times t 
characterize the rescaling intervals while r are the growth 
times within these intervals. 

Periods with LLLE< are associated with a return of 
skill in predictability experiments [llj ]. During these in- 
tervals the entropy production of the system is considered 
to be negative [ijj- 

IV. RESULTS 

The distribution of the LLLE is analysed for a stan- 
dard low resolution version and for variable vertical and 
horizontal resolutions. The validity of the Fluctuation 
Theorem is proven for a range of LLLE for a fixed growth 
time r in a low resolution experiment. 

In the standard low resolution version T21L5 of the 
model (total wave number 21 and 5 vertical levels) 7744 
degrees of freedom are present. Due to the numeri- 
cal efficiency this resolution is frequently used in long- 
term simulations since it produces a circulation with the 
characteristic properties of the observations (for exam- 
ple Hadley and Ferrel cells and mid-latitude storms [IB]). 
The LLLEs are determined for fixed growth time r = 100 
hours in all sampling intervals (tq = 10 days). Figure [T] 



shows the distribution of the LLLE values with a Gaus- 
sian fit (normalized to the maximum) . Heavy tails cannot 
be excluded but are unlikely. About 17% of the LLLEs 
are negative, hence the model is in a state of negative 
entropy production for a considerable amount of time. 
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FIG. 2. Quantiles (0.9, 0.99) of LLLE. The median (solid) 
corresponds to the global largest Lyapunov exponent 

To characterize the distribution for a wide range of 
growth times r the 0.9 and 0.99-quantiles are presented 
(Figure^]). The growth times r range from the minimum 
of 5 hours to the maximum given by the rescaling time 
To = 10 days. Thus Figure Q] corresponds to the r = 100 
hours slice in Figure [5J The median is the largest Lya- 
punov exponent which is by definition independent of the 
growth time r. For increasing r the quantiles approach 
the median indicating a narrowing distribution. For large 
r negative values become less likely and the probability 
of negative entropy production vanishes. 

The sensitivity of the LLLE distribution to the ver- 
tical model resolution is determined for fixed horizontal 
resolution T21 (Figure [3]) with numbers of levels ranging 
from 3, 5, 10, to 20. The number of degrees of freedom 
depends linearly on the number of levels: 4840, 7744, 
15004, 29524. Here Gaussian fits to the distributions are 
shown which are normalized by the maximum (compare 
Figure [T]) . All curves show a comparable mean but be- 
come broader for increasing numbers of levels. Beyond 
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FIG. 3. Gaussian fits of LLLE distribution in T21 with dif- 
ferent vertical levels L as indicated. 
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FIG. 4. Gaussian fits of LLLE-distribution for L = 5 levels 
and different horizontal resolution, as indicated. 



FIG. 5. Logarithmic probability ratio of positive to negative 
entropy production rates for increasing r = 10, 80, 160, 240h 
in T21L5 resolution. 



10 levels a limit is reached. 

For increasing horizontal resolutions and keeping the 
number of five levels fixed the LLLE distributions change 
differently (Figure [4} . The degrees of freedom for five 
levels are 4096 (T15), 7744 (T21), 16384 (T31) and 28224 
(T42). 

^From T15 to T21 the distribution shifts to positive 
means but becomes broader. Hence the broadening found 
here could be due to the increasing number of degrees of 
freedom as for the vertical resolution (Figure [3]) . 

/From T21 to T31 and T31 to T42 the main change 
is a shift to positive values while the variance increase 
is negligible (in contrast to Figure [3]). Therefore, the 
frequency of negative LLLEs decreases with increasing 
horizontal resolution. 

The increase of the horizontal resolution enhances the 
effective number of degrees of freedom by the simulation 
of smaller structures like mid-latitude vortices. In con- 
trast, the increase of the vertical resolution beyond 10 
levels does not yield additional degrees of freedom for 
the given horizontal T21 resolution. 

The validity of the FT is assessed with the T21L5 res- 
olution version. This FT relates the ratio of positive to 
negative entropy production rates £ with the entropy 
production and the growth time r 



P(-E) 



(2) 



Here P(£) is identified with the probability of observ- 
ing a LLLE of S. Figure [5] shows the ratio © for differ- 
ent growth times between r = lOh and r = 240h. These 
logarithmic ratios are clear linear functions of S. 



V. SUMMARY AND CONCLUSIONS 

The dynamic atmospheric model PUMA (Portable 
University Model of the Atmosphere) is used to assess 
the variability of predictability from a state space per- 
spective. The model has variable vertical and horizontal 



resolutions and is used here for L — 3 to 20 levels and 
I = 15 to 42 spherical wave numbers. The model de- 
scribes a standard dynamical core with diabatic heating 
and friction. In contrast to complex models, forcing and 
dissipation in PUMA are implemented as linear relax- 
ation processes. This keeps the system simple to concen- 
trate on the effects of the nonlinear dynamics. 

The local (time-dependent) largest Lyapunov exponent 
LLLE is estimated by the divergence of two close 200 
years trajectories, a reference and a perturbation; the lat- 
ter being rescaled in regular time intervals. The distances 
between both trajectories are measured by an Euclidean 
metric comprising all dynamic variables at all grid points. 

The main interest is in the frequency of negative LLLE 
which is associated with increased predictability (return 
of skill) and negative entropy production. Note that in 
contrast to the analysis of [8j who presented a regional 
analysis, the results in the present paper correspond to 
the global circulation. 

Furthermore the results arc of interest as they provide 
an example of a dynamical system which does not satisfy 
the conditions for a rigorous mathematical proof of the 
FT (it is not Anosov), but where it nevertheless seems to 
hold. 

For moderate resolutions below £ — 42 wave numbers 
the distribution of LLLE shows a considerable ratio of 
negative values (up to 20%) , hence negative entropy pro- 
duction and the violation of the second law. 

The increase of vertical and horizontal resolution shows 
different effects: The vertical resolution mainly increases 
the variance of the LLLE distribution while the horizon- 
tal resolution primarily increases the mean. These results 
suggest a method to determine an optimal combination of 
horizontal wave numbers and vertical levels in numerical 
models. The widening of the distributions for increasing 
vertical resolutions enhances fluctuations and can lead to 
an increased frequency of circulation extremes. 

The probabilities of the negative entropy production 
rates are in agreement with the Fluctuation Theorem [l|. 
These resolutions are still standard in many global warm- 
ing scenarios and even lower resolutions are still used for 
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long term integrations. The results show that numeri- 
cal models used in geophysical fluid dynamics may be 
considered as mesoscopic systems. 
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